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Wigner's' 1 ' suggestion to model the level statistics of complex quantum systems 
by random matrix models of very large dimensions is well known. In applications as 
diverse as two dimensional gravity' 2 ', topological field theory' 3 ' quantum chaos' 4 ' and 
mesoscopic physics' 5 ' different kinds of universal features emerge, all represented by 
particular types of random matrix ensembles. The analysis of these models seems to 
lead one to classical' 2 ' 3 ' and quantum integrable systems' 6 ' of very special mathematical 
structure and many unexpected relations make their appearance' 7 '. 

Random matrix models exhibit a large degree of robustness' 8 ' and the search for 
different kinds of universality has produced a wealth of results in two dimensional 
quantum gravity' 2 ' and in topological field theory' 3 ' . In quantum chaos and mesoscopic 
physics there are physical reasons to expect more universality classes than presently 
known' 9-12 '. This letter takes a small step in this direction. 

For reasons of concreteness and mathematical simplicity, we restrict our attention to 
a disordered system without time-reversal symmetry (e.g. in the presence of a magnetic 
field) and ask whether random matrix models can account for known departures from 
GUE (Gaussian Unitary Ensemble) level statistics' 10-13 '. On the one hand the class of 
matrix models with GUE level statistics has been shown recently to be very wide' 8 ' 14 ' 
and also to agree well with the properties in the energy region of extended states' 15 ', 
but, on the other hand, disordered systems have energy regions where the states are 
localized and the level statistics is Poissonian. To be sure, even localized states, once 
sufficiently close energetically, eventually repel each other because the exponentially 
small overlap can be outweighted by exceedingly small energy denominators. It is the 
typical scale of energy differences needed for level repulsion to be felt that changes 
dramatically in transition from extended to localized regimes. 

Recently, an attempt to incorporate both statistics in a random matrix model was 
made' 12 ': The basic idea was that in the localized regime the Hamiltonian prefers the 
site basis and hence the random matrix model was made to break the U (N) invariance 
which represents the lack of basis preference in the GUE ensemble. For sufficiently 
strong breaking, approximate Poisson statistics could be obtained. The explicit break- 
ing of U (N) invariance makes the mathematical treatment and the study of robustness 
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difficult. While the models we are discussing cannot be fully accurate representations 
of the mesoscopic physics since we will not introduce any explicit dependence on dimen- 
sion, we feel that exploring the possibility for constructing such an enlarged universality 
class is of more general interest and holds the potential for wider applicability than 
just to mesoscopic physics. 

Our main new observation is that basis preference, while needed, does not nec- 
essarily force the abolishment of U(N) symmetry. In a way vaguely reminiscent of 
spontaneous symmetry breaking, we argue that what is important is the preference 
for some basis but it is immaterial which basis this is. Averaging over all preferred 
bases we recover U (N) invariance and still manage to evade GUE statistics if a certain 
parameter is made strong enough. 

Imagine that the preferred basis is represented by an N x N unitary matrix V. This 
representation has N redundant phases and it is therefore convenient to consider instead 
the unitary matrix U = VDV^ where Dij = 5ije tdi , ( \6i\ < 7r, i — 1, 2, . . . N ). angles 
9i but For U to define uniquely a preferred basis we need the e j, s to be all different. 
In order to ensure this upon averaging we would like the weight to contain sufficient 
repulsion between the e l9j 's. Let the Hamiltonian be represented by a hermitian matrix 
H. For fixed U we consider the unnormalized distribution of H given by: 

Pu(H)d N2 H OC e -TrH> e - b Tr([U, H][ U,HV) d N* H (1) 

Here, d N2 H = Yl^dHaYl^j dRe(Hij)dIm(Hij) and b > 0. The 6-dependent term 
tries to align H with U and hence prefers the basis V. We now average over U with 
the invariant U(N) Haar measure dU . This induces enough repulsion between the 
eigenvalues of U to ensure that a unique preferred basis almost always exists. Our 
ensemble therefore is: 

P(H)d N2 H = C'(N, b) e~^ TrH2 [ J dUe 2bT ^ UHU ' H ^ }d N " H (2) 

C'(N, b) is a normalization constant. 
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The integral over U can be done' 17 ' and, going over to the marginal distribution for 
the eigenvalues of H , x iy we find: 

P(x 1 ,x 2 , ■ ■■x N )l[dx l = C(N,b)e-^ x * VeUj (e" 6 ^-^) 2 ) dx, 

i i 

C(N, tyVeUj (e-(P+h)tf+^)+*>xiXi} 7J dx . (3) 

i 

P(x 1 ,x 2 , ■ ■ -x N ) in Eq.(3) coincides' 18 ^ 20 ] with the diagonal element in the "coordinate 
representation (x)" of the density matrix of N non-interacting fermions of mass to, 
in one dimension, moving in a confining harmonic potential V(x) = \moj 2 x 2 at finite 
temperature, 

< Xi\exp{-(3{— H — )}\xj >J 

= M*f« a*fc(IM>)« + x|) + — J (4) 

where M is a normalization factor, flfiuj = arcosh(l + ^) and the length unit is fixed 
by rJ Y~ = y/l + 4b. This means that the mass and temperature are related by ^ = 
coth^. The GUE ensemble is obtained in the limit %uf3 — > 00, ^ — > 1 (implied 
by b — > 0, GUE statistics in Eq.(3) ) and complete disappearance of level repulsion is 
obtained in the limit hu(3 — > 0, ^ — > 00 (b — > 00, Poisson statistics in Eq. (3) ). 

When iV is large it is much more convenient to analyze the grand canonical version 
of the ensemble (3). Indexing the one particle states by a, denoting their energies by 
e a and the corresponding wave functions by ip a {x) one immediately can write down 
the expression for the local density n(x): 

- ( \ \- l^afo)! 2 

The chemical potential /i is fixed by AT via N = ^[e^ 6 " - ^ + I] -1 • 

For large N one can use semiclassical forms for ip a (x) as long as x is away from 
the edges of the support of n(x) and replace the sum over a by an integral. A simple 
computation gives then: 

n( x ) = (°° dp - (Q) 



Eq. (6) already exhibits both limiting situations, of GUE and of Poisson statistics: For 
b growing with N slower than N 2 , f3N — > oo, the gas becomes completely degenerate 
and one obtains the GUE semicircle law: 
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For b growing with N faster than N 2 , [3N — > 0, a gaussian distribution is obtained: 



-( \ N - 
v 7r 
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Correlations can also be easily calculated. Consider the 2-point connected correlation: 

x ^ x' 



5n(x)5(n)(x') = -\^2n a ip^(x)ip a (x') 



(9) 



with n r 



[e 



P(e a -n) 



+ l] -1 . Writing x = x + |Ax, x' = x - |Ax with < |Ax| << 



n(^) 1 we derive for (3N — > oo: 
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and for f3N -> 0: 



5n(x)5n(x') 



_ e -2b(Axf 
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The last two equations show that the model incorporates GUE and Poisson statistics. 

Consider now the intermediate case when, in the iV — > oo limit, j3N — > const = A. 
If the constant A is small, the average density of levels n(x) and the correlations are 
approximately the same as for f3N — > 0. The case A » 1 is more interesting: n(x) is a 



complicated function which, for small x, is proportional to \/A — x 2 and, for large x, is 
proportional to exp(—x 2 ) ("small" and "large" here mean, respectively, (A — x 2 ) » 1 
and —(A — x 2 ) >> 1). The correlation also changes when one moves from the center 
of the band towards the edges: at the edges (i.e. x » A) the levels are essentially 
uncorrelated, whereas at the band center (x ~ 0) the correlation function is given by 

(12) 



Sn(x)Sn(x') 



(n(x)y 



( SI 1 1 u ^ " 



iru 



ir 2 u/2A 



sinh(7r 2 M/2A) 



where u = n(x)Ax is the energy separation measured in units of the averaged level 
spacings. While a sequence of roughly A levels obeys GUE statistics, the correlations 
between more distant levels rapidly diminish. This is in sharp contrast with the case 
/3N — > oo (/3N — > 0) where, in the limit N — > oo, an arbitrary long sequence of levels 
obey the GUE (Poisson) statistics. In this sense the parameter j3N of our model can 
be identified with the dimensionless conductance g. g — > oo corresponds to a metallic 
behavior, g — > corresponds to an insulator and g — > const is identified with a mobility 
edge situation, where a third universal statistics, intermediate between the GUE and 
Poisson, is observed 10-12 . In that latter case AiVp xs N r /2A where AiVp describes 
fluctuations within an energy interval T with an average number of levels equal to iVr 
(Nr >> A is assumed). The linear relation between AiVp and N r is a consequence 
of the exponential decay of correlations in Eq. (12). In ref. 11 it was argued that 
a linear dependence is ruled out in the universal regime by a sum-rule expressing the 
conservation of the total number of levels, N. In our analysis of the matrix model we 
found that in the large N limit it is correct to pass from the canonical description - 
fixed number of levels or fermions - to the grand canonical one - arbitrary number of 
levels in the presence of a "chemical potential". Therefore, the overall level number 
conservation cannot affect universal features of the statistics of finite sequences of levels 
in the thermodynamic limit. 

A guess for the enlarged universality class of these models is: add more U(N) 
invariant function of U and H to the exponent in Eq. (1) and average over U. While 
strict U (N) invariance is useful, models with explicit breaking' 12 ^ are not excluded and 
can easily turn out to be controlled by U (iV)-symmetric fixed points. In this context 
it is amusing to note that if, in a reversal of roles, one integrates over H in (1) and 
views the resulting distribution of U as a generalized circular Dyson-type ensemble^ 16 ' 
one recovers a "circular" version of the eigenvalue distribution arrived at in Ref. [12]: 



Pe(9 1 ,...9 N )l[d6 t = C"(N,b)ll 
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1 +46sin 2 (^p) 
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For numerical purposes one may wish to stick to the canonical ensemble; expressions 
for n{x) and for ^"^^f ^ can be obtained in terms of infinite sums that converge fast 
as long as (3N does not grow to infinity t 20,21 !. For example, the exact formula for n(x), 
in Eq. (5), is: 

oo f oo N+r ) 

n(x) = \Mx)\ 2 S(-l) r 9 (JV " a)r+ ' r(r+1) II (! - rt" 1 (14) 

a=l \r=0 £=N+1 I 



where e a = hu(a -\),q = = [l + £ + ^ + ~ , and n^jv+i(l - (fl -1 = 1- 

Eq. (2) can be viewed as a degenerate case, in one space-time dimension, of "induced 
QCD" models' 22 '. Generalizations of the ensemble discussed here to the orthogonal 
or symplectic case are likely to produce technical complications similar to the ones 
encountered in the study of random unoriented two dimensional manifolds with or 
without matter[ 23 l 

In summary, the random matrix model proposed in eqs. (2,3) is an extension of the 
standard gaussian unitary ensemble with a level distribution that interpolates between 
Wigner-Dyson and Poisson statistics. U(N) invariance is preserved and the model is 
equivalent to a system of one dimensional fermions in a harmonic potential at finite 
temperature. Possible applications are to modeling of level statistics in quantum chaos, 
mesoscopic physics and disordered electronic systems. For example, in the latter case, 
the proposed ensemble can be driven from metallic to insulator behavior by varying 
the parameter b. 
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